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Abstract 



We consider a model of 2D quantum field theory on a disk, whose 
bulk dynamics is that of a two-component free massless Bose field 
X = (X,Y), and interaction occurs at the boundary, where the 
boundary values (Xb,Yb) are constrained to special curve - the 
"paperclip brane". The interaction breaks conformal invariance, 
but we argue that it preserves integrability. We propose exact 
expression for the disk partition function (and more general over- 
lap amplitudes (P\B) of the boundary state with all primary 
states) in terms of solutions of certain ordinary linear differential 
equations. 
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1 Introduction 



Lately, a class of 2D Quantum Field Theories (QFT) is attracting atten- 
tion from different points of view. These are the QFT which are Conformal 
Field Theories (CFT) in the bulk, but have non-conformal interactions at 
the boundary. In the simplest setting, the bulk CFT is placed inside the disk 
\z\ < R (where (z, z) are standard complex coordinates in the 2D space) and 
the non-conformal interaction occurs at the boundary \z\ = R. When the 
length scale changes from short to long, such theories typically interpolate 
between different conformal boundary conditions, and for that reason they 
are often referred to as the boundary Renormalization Group (RG) flows. 
Field theories of this type are of interest in dissipative quantum mechanics, 
where the bulk CFT theory plays the role P of the Caldeira-Leggett quan- 
tum thermostat [2]. They are also extensively used in studying the Kondo 
model and related problems c 3,'4j. Certain interest to such theories exists 
in the string theory 0E1IZ1- From the last point of view, of special interest 
are the "brane models", where the bulk CFT contains a collection of N 
free massless scalar fields X(z, z) whose boundary values X# = X||^ =R are 
constrained to some hypersurface 5] C M. N - the "brane" . Here and below 
we use the term "brane" loosely - the stringy on-shell brane must be a con- 
formal boundary condition, while we do not impose this requirement. The 
constraint gives rise to the boundary interaction which in general breaks the 
conformal symmetry: the shape of the "brane" changes with the renormal- 
ization scale. The RG flow equation can be computed perturbatively, in the 
limit when the curvature of the brane is small jSHH] . 

In this work we restrict attention to special model of this kind, which 
contains two scalars X = (X, Y), and the brane is a one-dimensional curve 
S C f 2 in the two-dimensional target space. We fix normalizations of these 
fields by writing down the bulk action 

A = - f d 2 z( d z xd- z x + d z Yd- z Y ) , (i) 

J\ Z \<R 

where d 2 z = dxdy. In this case, when S is one-dimensional, the RG flow 
equation was computed up to two loops in j^j . It turns out that this equation 
is satisfied by the following scale-dependent curve 

r cosh (^§) -cos (^-) =0, \Y b \<ttV^T2, (2) 

where n is real parameter independent of the RG energy scale E, but the 
coefficient r = r(E) "flows" with the scale E according to the equation 

K = nr n (1-r 2 ) , (3) 
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where k is inversely proportional to E: !<■=%-■ The proportionality coeffi- 
cient E* (the integration constant of the RG flow equation) sets the "physical 
scale" for the model: physical quantities, like the overlap amplitudes © be- 
low, will essentially depend on the dimensionless combination E*R. In what 
follows we always take the normalization scale E equal to R , so that k in 
the left-hand side of (j3J) coincides with this combination, 

K = E*R . (4) 

Eq.© then relates r to the radius R. The curve (J2J) is sketched in Fig.l. 
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Figure 1: The paperclip shape for r = 2 x 10 . 

It has a paperclip shape, and we will refer to this theory as the "paperclip- 
brane model", or simply the "paperclip model". At short scales, i.e. at 
R — > 0, the coefficient r becomes small, and the paperclip grows long in 
the horizontal direction, with its length ~ -7= log (-), while its width ap- 
proaches the constant tt \fn + 2. 

Strictly speaking, the two-loop paperclip solution ©,© applies only 
when n>l, for otherwise the curvature at the round ends of the paperclip in 
Fig.l is not small, and higher loops and perhaps nonperturbative corrections 
should not be neglected. For similar reason, Eq.(J2J) can be taken literally 
only at sufficiently small length scales, where r is well below 1. For larger 
length scales a nonperturbative description of the paperclip-brane model is 
needed. At the same time, there are reasons to believe that existence of con- 
sistent local boundary conditions is not limited to the perturbative domain 
u> 1, and in this sense the "paperclip" model can be given appropriate 
nonperturbative definition for all n > 0. It is the problem we address in this 
work. 

Nonperturbative description of a local boundary condition can be given 
in terms of the associated boundary state ^01- The boundary state | B) is 
a special vector in the space of states Ti of radial quantization of the bulk 
CFT which incorporates all the effect of the boundary at \z\ = R. In our 
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case H is the usual space of states of two-component uncompactified scalar, 

H = j J=p®?p, (5) 

where J-^> is the Fock space of two-component right-moving boson with the 
zero-mode momentum P = (P, Q). The overlap (v\B) with any | v ) E T~L 
coincides, up to the factor j{ c /6- 2 ^ v ^ w jth the (unnormalized) expectation 
value ( O«(0, 0) )disk of the field O v (z,z) associated with \v) through the 
standard operator-state correspondence of the bulk CFT. In particular, the 
overlap ( P | B ) with the Fock vacuum | P ) relates to the unnormalized one- 
point function of associated primary field inserted at the center of the disk, 

(j™(0,0)) disk = R^- F2 /*(P\B), (6) 

and i? 1 / 3 ( | B ) is the disk partition function 1 . When the boundary con- 
dition is not a conformal one, the overlaps depend on the combination (0J) ; 
to stress this fact we will use the notation 

(P|fl) = Z(P|«). (7) 

In this paper we propose exact expression for the function Z(P\k) in 
terms of solutions of certain ordinary linear differential equation. Our pro- 
posal is summarized in Eas. ()9U|) - (|95|) in Section \7. II below. The motivations 
and evidence for this proposal are presented in turn in the main part of 
the paper. Here we would like to outline the logical structure of our argu- 
ments, and their relation to known exact results in other models of boundary 
interaction. 

First, we are going to argue that the paperclip-brane model is integrable. 
More precisely, we will argue that this model fulfills the paradigm of bound- 
ary integrability as described in Ref . ^2] • We will show that the holomorphic 
sector of the bulk CFT (^Q) contains special set of commuting local integrals 
of motion { I s ; s = 1, 3, 5 . . . }, and that the paperclip boundary constraint 
preserves the combinations I s — I s , where I s are the antiholomorphic coun- 
terparts of I s . In the radial quantization picture, where I s are hermitian 
operators acting in the space J-p, such that 

[I, , I S ' ] = (8) 

x We assume that the exponential field in @ is defined according to the usual CFT 
conventions | 11| . The same will apply below to all composite fields in the bulk of the 
disk. Namely, unless specified otherwise, all bulk composite fields appearing below are 
understood as the normal-ordered products, taken with respect to the standard Wick 

paring z)X"(w, w) = —5^" log \z - w\. 
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for all s, s' , the precise form of the last statement is that the paperclip-brane 
boundary state | B ) satisfies the equations 

(I s -I s ) \B} = (9) 

for all s = 1,3,5... . The first integral Ii essentially coincides with 
where Lo is the zero- mode Virasoro generator, and so the first of the equa- 
tions is simply the statement of the translational invariance of the boundary 
condition (j2j). Few further operators of the "paperclip series" I3, I5 . . . are 
described in Section |S] below. General explanations about the local meaning 
of the equations @ can be found in Ref. [P2] . 

The equations © imply that the boundary state can be written in terms 
of the simultaneous eigenstates | a , P ) € Tp of the operators I s as follows, 

\ B )= I J2 B ^ P ) |a,P)®h7^), (10) 

(where the over lined states are the corresponding eigenstates of { I s } in Tp) 
with some c- number coefficients B a (P). Also, one can observe (see Section |2] 
below) that the expressions for the operators I s do not involve any particular 
scale, hence the eigenvectors | a , P ) do not depend on the scale E*, while all 
the dependence on the parameter (jlj comes through the coefficients B a (P) 
in (|l(Jj) : B a (P) = B a (P \ k). This structure is typical for all integrable 
models of boundary interaction; it suggests interesting general relation ^H] 
between the integrable boundary states and commuting operator families a 
la Baxter |14j . Assuming that the integrals of motion I s are "resolving", i.e. 
that for generic P all eigenspaces of { I s } in T-p are one-dimensional (which 
seems to be the case), one can associate to the state (fTU|) an operator B(k) 
( the "boundary-state operator") acting in the space J p Tp, 

»(«)= / y>a(P|«) |a,P)(P, a\. (11) 

In writing (|11|) we have also assumed that the eigenstates | a , P ) are chosen 
orthonormal, 

(P,a|a , ,P , ) = ^ ^(P-P'). (12) 

Obviously, the operators M(k) with different values of the "spectral param- 
eter" k commute, 

[B(/s),B(«')] =0, (13) 
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in suggestive similarity with Baxter's operators. And indeed, it was observed 
in ^3] that the appropriately defined Baxter's operators of minimal CFT 
can be identified with the boundary-state operators associated with certain 
integrable boundary conditions, and this relation seems to be of a general 
nature [To1lltjlll7j . Two more observations about Baxter's operators of CFT 
will be important for our arguments. 

First, not only the operator B(k) commutes with the local integrals {I s }, 
but in some cases it admits useful asymptotic k — > oo expansion in terms 
of these integrals. Assuming that the index s labeling the local integrals I s 
coincides with their scale dimensions, this expansion can be written as 



where Mm is the boundary state operator associated with the infrared fixed 
point, and E* relates to k as in (|4*)l. The dimensionless constants C s de- 
pend only on accepted normalizations of the integrals of motion I s . In all 
known examples, the asymptotic expansion (|14|) is valid when the bound- 
ary flow associated with B(k) ends in the "trivial" infrared fixed point, the 
one which in Cardy's classification ^H] corresponds to the identity primary 
state; otherwise, the k — > oo expansion of logB(K) involves also non-local 
integrals of motion of fractional scale dimensions |15j . Conversely, when 
an integrable boundary interaction flows to the "trivial" fixed point, the 
structure (|14jl is natural to expect, because the "trivial" conformal bound- 
ary condition features no local boundary fields but the descendents of the 
identity operator jTH]. 

Another important observation belongs to Dorey and Tateo, who have 
shown that, in the case of minimal CFT, the vacuum- vacuum matrix ele- 
ments of Baxter's operators (i.e. the overlap amplitudes analogous to Q) 
can be written down explicitly, in terms of solutions of certain ordinary lin- 
ear differential equation |19j (see also HO]). Albeit still rather mysterious, 
this relation has proven to be rather general and very useful. The differential 
equations have been identified for quite a few different classes of Baxter's 
operators, where the relation can be either proven |20|I16| or verified against 
known "TBA systems" of functional relations. In our case, for the paperclip 
model, the associated TBA system (which is related to the TBA system of 
the "sausage" sigma model is too complex to take it as the starting 
point. On the other hand, a wealth of data about the paperclip boundary 
state can be obtained through perturbative analysis in the weak coupling 
domain, and by looking into various limiting cases of the model; Sections 
01 through HO below are devoted to these tasks. Then one can simply try to 




(14) 
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figure out the differential equation which would reproduce all these data. 
This is the approach we have taken in this work. In Section [7] we present 
special second-order differential equation, Ea. ()9U|) . and propose the relation 
(|95j) expressing the amplitude (Q) in terms of certain solutions of that dif- 
ferential equation. Then we show that this expression reproduces expected 
properties of the paperclip model in all fine details. 

The last thing we would like to mention before entering the details is 
the circular limit of the paperclip. In the limit when n — > oo and r is taken 
sufficiently close to 1, so that g -1 = n (1 — r 2 ) remains finite, the paperclip 
curve © becomes a circle, 

X| = i , with k = g~ l , (15) 

and we refer to this limiting case as the "circular brane" model. This model 
is important in view of its applications in Condensed Matter Physics |22| 
I23H241I55] . in particular in relation to the problem of Coulomb charging in 
quantum dots. The results for this special case were previously reported in 
Ref. [HI- 

2 Elementary Paperclip 

In this section we list few elementary properties of the paperclip model, 
mostly to prepare the notations for the future discussion, and to set the 
stage altogether. 

2.1 Functional integral 

The one-point function © can be represented in terms of the functional 
integral as follows, 

( e ip ' x (0, 0) > disk = J VXVY e [px+ ' lQY (0, 0) e~ A ^ , (16) 

where P = (P, Q), and the integration variables X(z, z), Y(z, z) are assumed 
to obey the constraint (J2J) at the boundary \z\ = R. In view of the compact 
nature of the paperclip ((2J) , the integrand in (|16|) is bounded for any complex 
P and Q, and therefore the overlap (P | B ) in Eq.® is expected to be an 
entire function of these variables. For this reason it is useful to think of the 
momenta P, Q as the complex variables. For instance, when P and Q are 
pure imaginary, some insight can be gained by making a shift of integration 
variables 

X^X + iPlog^, (17) 
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in the functional integral (|16|). which brings it to the form 

(e ip ' x (0,0)) disk = ir p2 / 2 J OTWe-- 4 ^'^^^, (18) 

where the parameters (P, Q) play the role of external fields coupled to the 
boundary values Xb = (Xb,Yb) via the boundary action 

A B = -l P- {PXb + QY B ) (z) . (19) 

J\z\=R 27TZ 

2.2 Topological sectors and instantons 

The functional integral (jlfij) is taken over all field configurations X(z, z) 
inside the disk satisfying the paperclip constraint © at the boundary. Since 
topologically the paperclip (J2J is a circle, the configuration space splits into 
the topological sectors characterized by the integer-valued winding number 
w £ Z. This is the number of times the boundary value winds around 
the paperclip when one goes around the circle |z| = R. As usual, existence 
of the topological sectors allows one to generalize the paperclip model by 
introducing the topological #-angle, so that the contributions from different 
sectors are weighted with the phase factor e l9w . The boundary state | B) 
would then acquire the ^-dependence. In this work we do not attempt to 
study the full ^-dependent paperclip model, restricting attention to the case 
9 = 0. Nonetheless it is important that even in this case the partition 
function Z( P | k ) can be written as a sum 

Z(P|«) = X;Z< tt >(P|«), (20) 

of the contributions from all the winding sectors. 

The decomposition Q2U|) is most useful in the semiclassical domain, where 
the contributions from different topological sectors are easier to sort out. 
The semiclassical approximation is valid when the overall size of the paper- 
clip (J2J) is sufficiently large. This is the case when n 3> 1 and log (^) is 
not too small, so that r"Cl Note that according to © one has to have 
sufficiently small R in order to meet the last condition; in other words the 
semiclassical domain corresponds to large n and sufficiently small length 
scales. In the leading semiclassical approximation the contribution from 
given topological sector is dominated by the classical solutions minimizing 
the action (^Q) in that sector - the instantons. It is not difficult to find the 
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instanton solutions. At large n one should not distinguish between n and 
n + 2 in (|2"|). and the action (^Q) can be written as 

Aa "' = ^L R Az — ^ — ' <21) 

where the complex- valued fields, 

= e (^+i^)/v^ j £/* = e (X-iY)/Vn (22) 

satisfy the paperclip constraint which in this context is best written in the 
form 

[rU B - 1) {rU* B - 1) = 1 - r 2 , (23) 

with Ub = U\\ z \ = r. In a given topological sector w the action (|21|) is 
bounded from below, Classical > M A/2ir, where A = 2tt log (^) is the 
area of the paperclip (J2J) . This bound is established by standard manipula- 
tions 27 j, with the winding number u; written in the form 



d,Ud*U* - drUdM* 



w = y I - - , , - - ■ (-'» 

For u; > one finds 



A J ]Z]<R UU* 



( w ) _ wA _ n 
class 2vr ~ 7T 



| d 2 z 
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2 >0; (25) 



this equation also shows that the minimum is achieved when U = U(z) is a 
holomorphic function of z inside the disk, satisfying the constraint (|23j) at 
the boundary. Since Eq.(J221) defines a circle, such functions exist and have 
the form 

where are arbitrary complex numbers inside the disk, \a,k\ < R, and 4> is 
an arbitrary phase 2 . One easily checks that U(z) defined by this equation 
is indeed holomorphic for \z\ < R as long as r < 1. Similar expression, with 
z replaced by z, can be obtained in the case of negative w. The parameters 
a,i and (j) are the moduli of the u>-instanton solution. Roughly speaking, the 

2 Eg. 1261 generalizes the instanton solutions of the "circular-brane" model, which were 
found in 
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deviations R — |oj| represent sizes of the individual instantons, while the 
phases of a{ correspond to their locations along the boundary circle \z\ = R. 
For any tf-instanton configuration 

ex P (-4ri)=^ nH > (27) 

i.e. the contributions from the sectors with are suppressed by small 

factors r n l ro l 

Nonetheless, instantons in the paperclip model give rise to a subtle effect 
generally known as the "small-instanton divergence" (see |3E,)- Evaluation 
of the instanton contributions involves integration over the moduli space, 
and this integral diverges logarithmically when the parameters a{ approach 
the circle |aj| = R. The divergence is clearly seen in the one- instanton 
calculation, see Section Hi 21 b elow . This divergence can not be absorbed into 
the renormalization of the parameters of the curve © ; one has to introduce 
an unusual boundary counterterm proportional to r n . Similar phenomenon 
is known to exist in the 0(3) sigma model {301131] and in circular case of 
our model (see e.g. |26| and references therein). One important effect of the 
small-instanton divergence is that it gives rise to a weakly singular factor in 
the partition function 0, 

Z(P\k) = (Ck) 2k Z(P\k), (28) 

where Z admits normal perturbative (in the powers of r) and "instanton" 
expansion (in the powers of k). The factor C in (|28j) represents ambiguity 
which derives from the small instantons - it depends on specific regulariza- 
tion of the small-instanton divergence. We will say more about this effect 
in Section 13.21 below. Here we just observe that the ambiguity affects the 
extensive energy of the paperclip boundary, i.e. the coefficient / in the 
asymptotic of Z at large R, 

logZ(P|«) /(P)k as k^oo. (29) 

In what follows we assume that the ambiguity is fixed by imposing the 
normalization condition 

/(0) = 0, (30) 

which is always possible to achieve by appropriate adjustment of the constant 
C in QS) . 
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3 Semiclassical Paperclip 



Here we evaluate leading semiclassical contributions to the overlap © by 
direct calculation of the functional integral ()16|) in the saddle-point approx- 
imation. This will give some intuition about its structure. We will consider 
contributions from the zero- and the one-instanton sectors only. 

3.1 Zero winding 

Let us first assume that the parameters P, Q are sufficiently small so that 
vertex insertion in is "light", i.e. it has no appreciable effect on the 
saddle-point configurations. We write 

(P,Q) = ^ (p,q) (31) 

and assume that p, q ~ 1. Then, at zero winding the action is minimized 
by the trivial classical solutions X(z,z) = Xo, where Xo = (Xq,Yo) is an 
arbitrary point on the paperclip (J2J). Therefore 

4a S = s 0) = / dX(x ) e np^+ q y )/v^ t (32) 

J paperclip 

where the integration measure dA^(Xo) is determined by integrating out 
fluctuations around the classical solution in the Gaussian approximation. 
Of course there is no need of actually evaluating this Gaussian functional 
integral to figure out the answer. If one writes X(z, z) = Xo + 5~X.(z, z), and 
splits the fluctuational part <5X into the components normal and tangent to 
the paperclip at the point Xo, the components are to satisfy the Dirichlet 
and the Neumann boundary conditions, respectively. Therefore one just has 
to take the product of the known (see e.g. Appendix to jH2]) disk partition 
functions with these two boundary conditions. As the result, dA^(Xo) = 

f|d^(X ), where d^(X ) = ^{dX ) 2 + (dY ) 2 is the length measure of 

the paperclip, and go = is the "boundary degeneracy" jSS] associated 
with the Dirichlet conformal boundary condition 3 . Certainly, this argument 
ignores the curvature of the brane at Xo, but its proper account leads to 
the one- loop part of the renormalization of the parameter r, as described by 
Eq.@ (see 0). 

3 The definition is as follows: gn — (P\ Bd ), where | Bd ) is the boundary state of 
uncompactified boson X with the Dirichlet boundary condition Xb = 0, and the primary 
states \P) are delta- normalized, (P \ P') = S(P - P'). 
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Of course there are many ways to evaluate the integral ()32j) . but the 
instructive one is to change to the variable 

Uo = e (x 0+ iY )/VE ; u * = V<£L i (33) 
where the expression for Uq follows from H23[); this brings (|32[) to the form 

4r s 0) = AVW^j* (34) 



2vri 



r 



The integration here is taken over the circle |r-J7o — 1| = y/l — r 2 (see Eq.(|23|l) 
which goes around the branching points Uq = r and Uq = ^. By closing the 
contour on the branch cut between these two points the integral is brought 
to a standard representation of the hyper geometric function, 

4as = s 0) = 9 2 d Vn(l-r 2 ) r-^ 2 F X (A + q - ip, \ - q - ip, 1; 1 - r 2 ) . (35) 

It is also instructive to rewrite this result in another form. The integration 
contour can be represented as the combination of two contours, both going 
from ^0 = back to Uq = 0, one around the point Uq = r, and another 
around Uq = -. Again the hyper geometric integrals emerge, and one obtains 

Z cia^ 0) = B clas S (p, q) FdassiP, q\r) + B class (-p, q) F dass (-p, q I r) , (36) 
where 

r(±-q-ip)T(± + q-ip) 

and 



Fcto*(p, q\r) = Vl-r 2 2 F l (±+q + ip,±-q + ip, 1 + 2ip; r 2 ) . (38) 

In this form the nature of singular behavior at r — > is more explicit. Let 
us make a (trivial) observation that the poles of -B c i aS s(p> q) i n the variable p 
in the first term in (|36|) are cancelled by the poles in the higher terms of the 
hypergeometric series -F c i a ss i n the second term, and vice versa, in agreement 
with the statement that the full partition function is an entire function of 
P and Q. 

Although the above result was derived under the assumption that P, Q 
are small, it needs little fixing to become valid for much larger values of these 
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parameters. When (P, Q) become as large as y/n the vertex insertion in (|16|) 
is "heavy" , i.e. it must be treated as a part of the action, and it affects the 
saddle-point analysis. In this case the form (|18|) of the functional integral is 
more convenient. The saddle-point configuration(s) is still a constant field, 
X(z, z) = Xo (this time we are talking about the shifted field, Eq.(|17jl). but 
now Xo is not an arbitrary point on the curve @, but has to extremize the 
boundary action (|19j) . 

•A s [Xo] = -iP-X • (39) 

There are two solutions of the saddle-point equation. They are easily visual- 
ized in the case of pure imaginary P and Q, when the associated saddle-point 
field X is real valued: the points Xo are located at two opposite sides of the 
paperclip (|2*|). one corresponding to the minimum, and another to the maxi- 
mum of the action (|3tJ|) . Of course the minimum dominates, but it is useful 
to keep both in mind. The saddle-point action plus the Gaussian integral 
over the constant mode produce nothing else but the p, q — ► oo asymptotic 
form of the expression (|35j) - after all, this asymptotic of the integral (|34|) 
is controlled by the same saddle points. One can observe that if one splits 
the constant-mode integration into two parts, as was suggested in deriving 
(|36j) . the parts receive contributions from different saddle points - one from 
the "minimum" and one from the "maximum" one; this is why one of the 
terms in (|36[) becomes negligibly small at QmP ~ y/n. What makes the 
difference at (P, Q) ~ y/n is the proper treatment of non-constant modes. 
One writes 

X(z, z) = X + t 5X t (z, z) + n 5X n (z, z) , (40) 

where Xo is the position of the saddle point on the paperclip, and to and 
no are unit vectors tangent and normal to the paperclip at this point. Then 
for small 6Xt the boundary constraint (j2J reads 

SXn = -^Ttt=^y cosh 6X ? + °( 6X ^ ' ( 41 ) 

2yn(l — r z ) v 
and, up to higher-order terms, the boundary action l)19j) can be written as 

P f dz 

AB-MXctezj-fj^tX}. (42) 

with the coefficient P r proportional to the normal component of the "exter- 
nal field" P times the curvature in (|41[). Explicitly, 



Pr ^ . (43. 
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where we assume that the branch of the square root is chosen in such a way 
that QmPr > 0; then the sign plus (minus) in (|42[) applies to the "minimum" 
("maximum") saddle point. Thus, while to the leading approximation the 
normal component 5X n still can be treated with the Dirichlet boundary 
condition, the "boundary mass" term in (|42|) has to be taken to account 
in evaluating the contribution from 5Xf. Note that for ~ 1 (which we 

assume) the energy scale associated with this "boundary mass" is ~ R~ , 
so that the use of the renormalized paperclip parameter r defined as in © 
is still appropriate. Using the well known boundary amplitude of the free 
field with quadratic boundary interaction (|42|) 33j, one finds that Eg . (|36|) 
would apply to the case of (P, Q) ~ y/n as well if one puts corresponding 
additional factors in the two terms in (jHEJ), i-e. replaces -B c iass (p, q) there by 

B dass (P, Q) = B dass (^ P,f Q)r(l-|); (44) 

of course in this case one can use the p, q — > oo asymptotic forms of the 
factors (HSl) and (gHJ) . 

3.2 One-instanton sector 

Here we consider only the case of the "light insertion", i.e. assume that P 
and Q are of the order of as in (j31|) . Then the saddle point configurations 
are the instanton solutions (|26j) . All the instanton solutions with w = 1 are 
generated by SU(1, 1) maps from the "canonical" solution Uq{z): 

with the moduli o, a* and <p appearing as the parameters of the map. The 
saddle-point contribution from the sector w = 1 to the functional integral 
(|16|) is the integral over the moduli space, 

4ll 1} =r n f dM(a, a*, 0) ([/ ) ip+<? (U^" , (46) 

where the factor r n is due to the saddle-point action ()27|) and C/o stands 
for the solution (j45[) evaluated at the center of the disk, C/o = U(0) = 
Uq(— e^a). The integration measure dA1(a, a*,^) should be computed as 
usual by integrating out the Gaussian fluctuations around the instanton 
solution. The calculations are analogous to those in [251 but tedious, we will 
present them elsewhere. The resulting form of the expression 1)46(1 is 

n f d 2 a 9 2 D {n(l-r^ ^ ^ i +ip+q i +ip _ q 
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where d 2 a = d(5Rea) d(3mo), and again Uq stands for U(0) = Uq(— e 1 ^ a). 

The small-instanton problem is explicit in (|47|) , The integral over \a\ 
diverges when \a\ — ► R. This limit corresponds to small instantons - when 
\a\ is close to R the function U(z) in (|4*5|) is almost constant (equal to 
the "background" value Uo = U(0)) everywhere in the disk except for a 
small neighborhood of the point z = a, where the instanton is localized. 
The integral ()47|) requires regularization. Any regularization is essentially 
an instruction to suppress contributions from the instantons of the sizes 
A < A , where A is the ultraviolet cutoff of the theory. The integral (|47|) 
thus produces a A-dependent term, 

= « log M4:: 0) + finite. (48) 

Similar divergences appear in the multi-instanton sectors, and by the usual 
arguments (which also take into account the "instanton - anti-instanton" 
configurations) the divergent contributions exponentiate, leading to the over- 
all A-dependent factor 

exp (2k log(A-R)) (49) 

in the partition function Z(P,Q \ k). The term 2k log A in this exponential 
can be absorbed by local boundary counterterm 

2i E*R I i^log(cl-), (50) 

J\z\=R l7:z 

leaving behind the singular factor shown in (|28|) . Note that the counterterm 
(|5Uj) has singular dependence on the paperclip parameter r - that is what 
makes it distinct from usual perturbative counterterms. 

All possible regularizations of the integral (|47|) essentially differ in the 
notion of the "size" A of a small instanton. Here we assume a natural 
definition 

R 2 - aa* 

A= ^mr- <51) 

which respects all evident symmetries of the paperclip © 4 . With this choice, 

4 Of course this choice is not unique. Generally speaking, any choice different by a 
factor, A — » C(Xo, Yo) A, where Xo, Yo is a "background point" related to Co as in 13311 . 
is acceptable as long as C(Xo, Yo) respects the symmetries. Obviously, this ambiguity is 
related to the possibility to add counterterm 15UH with C being a local function of the 
boundary fields, C = C(Xb,Yb)- Therefore in our case (unlike the circular-brane model 
where the symmetry restricts the ambiguity to a constant C) regularization of the small- 
instanton divergence does have certain effect on the physical content of the theory. From 
this point of view 151H is rather a part of our definition of the paperclip model. We will 
see that this choice is consistent with the integrability of the model. 
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the easiest way to handle the integral (|47|) is to insert regularizing factor A 2e 
into the integrand. Then, after subtracting - [CE if )~ e Z^~^' (which is the 
job of the counterterm (|5U|)). the limit e — > can be taken. The integration 
can be handled in a way similar to that in the previous subsection, and the 
result for the finite part in (|48|) is 

\gl (n(l-r 2 ))i r"- 2i "x (52) 
^ 2F 1 (l+q-ip + e,±-q-ip + e,l + 2e;l-r 2 ))\ t=0 . 

And similarly to (|36|) . this part of the one-instanton contribution can be 
written as the sum of two terms behaving at small r as r nziz21p times a power 
series in r 2 . 

4 Ultraviolet Paperclip 

As was mentioned in Introduction, at short length scales the parameter r 
in |2J) becomes small and the paperclip grows long in the X-direction. In 
the limit r — * it can be regarded as a composition of two "hairpins" , as 
depicted in Fig. 2. 



Figure 2: The paperclip formed by a junction of two hairpins. 

Qualitatively, one expects that if QmP is not too small, the R — > 
behavior of the overlap (P\B) is controlled by one of the hairpins, at the 
left or at the right in Fig. 2, depending on the sign of QmP, with some 
crossover at small QmP. It is useful therefore to describe the boundary 
state of the hairpin brane in some details. 

4.1 The hairpin 

For definiteness, we will discuss here the hairpin in the left-hand side of Fig. 2; 
we call it the left hairpin. The right hairpin is the X — > —X reflection of 
the left one. The (left) hairpin is described by the equation 

^~ ex P(-S)- cos (^f)=°- ( 53 ) 
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The hairpin is a brane in usual stringy sense - the boundary condition (|53jl 
is a conformal one. More precisely, the curve Q53JI satisfies the RG flow 
equation with 

f (54) 

where k = E*R. We have checked this up to two loops, but it likely holds 
to all loops 5 . The equation Q54|) shows that when the RG "time" t = log k 
increases the hairpin "flows" uniformly to the right with no change of shape. 
It can be made into an RG fixed point in a precise sense of the word by an 
appropriate redefinition of the RG transformation, namely by supplementing 
it with a simple field redefinition (X,Y) — > (X + -^=,Y). This corresponds 
to introducing a modified energy-momentum tensor 

t = -d z xa z x-d z Yd z Y + -^d 2 z x, (55) 
f = -d z xd- z x-d- z Yd z Y + ^d 2 z x, 

or, in a stringy speak, to introducing a linear dilaton -D(X) = -j= X (for 
the problem at hand such modification, in whatever speak, has no effect on 
the physical content of the theory). With this, the boundary state \B) C 
associated with the hairpin brane enjoys the conformal invariance in the 
usual form, 

[z 2 T(z)-z 2 f(z)} lzl=R \B) c =0. (56) 

Moreover, the hairpin brane has an extended conformal symmetry generated 
by higher spin currents - the VF-symmetry. The generating currents W s {z) 
have spins s = 2,3, 4. . . (with W2(z) = T(z)), and can be characterized by 
the condition that they commute with two "screening operators", i.e., 

j dw W s {z)V±{w) = 0, (57) 

where 

V±=e V^X±iV^+2Y ? (58) 

and the integration is over a small contour around the point z; vanishing 
of the integral (|57|) implies that the singular part of the operator product 

As usual, the higher loop corrections are scheme dependent. The statement is that a 
scheme exists in which Ea. l|53^ is valid to all orders in the loop expansion. 
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expansion of W s (z) V±(w) is a total derivative d w (. . .). This condition fixes 
W s (z) uniquely up to normalization and adding derivatives and composites 
of the lower-spin IV-currents. For instance, the first holomorphic current 
beyond can be written as 

W 3 = ^±l{d z Yf + n{d z X) 2 d z Y + 

ri\/n o (n + 2)\/n „ „ 9 n + 2 o , 

where the ambiguity in adding a term proportional d z T is fixed by demand- 
ing that (|59[) is a conformal primary. One can notice that the current W3 is 
antisymmetric under reflection 1" — ► — Y", but has no symmetry with respect 
to X — ► —X. The higher currents W4, W5, . . . can be found either by a 
direct computation of the operator product expansions with the screening 
exponentials (|57j). or recursively, by studying the singular parts of the op- 
erator product expansions of the lower currents, starting with W^{z)W 3 {w) 
and continuing upward. Thus, the product W${z)Wz{w) contains singular 
term ~ [z — w)~ 2 which involves, besides the derivatives d z T and the com- 
posite operator T 2 , the new current W4. Further operator products with 
W4 define higher W's, etc. In this sense the VF-algebra is generated by 
the two basic currents T and W3. This M^-algebra itself was previously 
described in Ref. 3"> . We also mention in passing that when n is a posi- 
tive integer the VF-algebra thus defined closes within n — 1 lowest currents 
W2, W3, . . . W n and coincides with the VFA n _i-algebra [36|l37| with the cen- 
tral charge c(n) = 2 + — . 

The hairpin boundary condition is invariant with respect to this TV- 
algebra, i.e. 

[z s W s (z)-z s W s (z)] lz]=R \B) c = 0, (60) 

for all the W-currents. To verify this statement by direct computation with 
the constraint Q53|) is not that easy. We have checked that it holds true 
in the classical limit n — > 00, where the calculation reduces to verification 
that the difference z 3 W3 — z 3 W3 vanishes at the boundary in virtue of the 
classical equations of motion of the hairpin model. In fact, in the classical 
case this condition turns out to be very defining, in the sense that it fixes 
the shape of the hairpin uniquely, up to overall scale of the fields (X, Y) 
(in the classical case n — > 00 one should not distinguish between yfn and 
y/n + 2 in ) . It is plausible that Eq. (jB0|> can be verified order by order in 
the loop expansion, but we did not perform any part of this task. Instead, 
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we have taken an attitude that Ea. (|6U|) . combined with Cardy's consistency 
condition JH! and qualitative properties of the hairpin lj53j). can be taken 
as the definition of the quantum hairpin 6 . 

In view of (|6U|) the hairpin boundary state can be written as the combi- 
nation of the Ishibashi states | Ip ) associated with the above W-algebra [HH] . 
An individual Ishibashi state 

|/p>= [1 + V( p ) X^X^ + ...] |P> (61) 

is the solution of the equations (|6U|) in the Fock space T-p <8> T-p with 
the zero-mode momentum P = (P,Q). Here X_, = (X_&, Y1&) are the 
bosonic creation operators associated with the oscillatory modes of the fields 
X^ = (X,Y). The solution is unique, i.e. the equations (|6U|) determine 
uniquely the coefficients A^iP), as well as all the higher-level coefficients; 
for instance, the amplitudes explicitly written in (|61j) are (we rather 
write down certain linear combinations which are shorter and more sugges- 
tive) 

2i (n + 1) . , 

Axx - A YY ± - \ \ A XY = 62 
y/n(n + 2) 

1 + n 1 + 2i^nP=F 2 v / rT+2Q 



l + n-2\^RP 1 - 2i^ P T 2 v / ^+~2 Q 



Axx+Ayy = . —^- x ( 63 ) 

1 + 4n 2 P 2 + 4(n + 2) 2 Q 2 
(1 - 2i^nP - 2V^+~2 Q)(l - 2i^-P + 2 v / ^+f<2) ' 

The hairpin boundary state is a superposition of these Ishibashi states, 

\B) C = J dPdQ B(P,Q) |J P ), (64) 

where B(P,Q) coincides with the vacuum overlap of the hairpin boundary 
state, B(P,Q) = (~P\B) C . Its exact form can be figured out from the 

6 In this approach one still has to check that this definition agrees with the perturbative 
definition of the hairpin brane through Eg. 15311 . We have checked that Eqs. l62t and 1651 
below are in agreement with the one-loop calculations on the disk with the boundary 
condition 1531 . 
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anticipated singularity structure, 

B(P,Q)=g 2 D (£)^x (65) 

v^r(-iy^p)r(i-iP/v^) 

r(i_^+2f _i^f)r(i + x ^+2f -Vnf) ' 

and then the overall expressions (|64[) . (|65j) can be checked against the con- 
tinuous version of Cardy's consistency condition ^B]- Here we skip this 
part in view of the recent paper which seems to contain closely related 
calculation for the case of Dl-brane in the "cigar" sigma- model. 

Note that the ^-symmetry ()6U[) of the hairpin model, together with its 
defining relations , strongly suggests the possibility of "dual" description 
of the hairpin model, similar to the duality between the conformal cigar 
sigma- model and the "sine-Liouville" model [101 1^ • Indeed, such dual 
description is possible, but we did not yet elaborate all its details; we hope 
to come back to this issue elsewhere. 

As was mentioned above, the right hairpin is the X — > —X reflection of 
the left hairpin described above. In particular, for the boundary state \B) D 
of the right hairpin, the overlap amplitude ( P | B } D equals B(—P, Q), and 
the whole boundary state of the right hairpin satisfies the Vl^-symmetry 
conditions 

[z s Wp\z) -z s Wp\z)} lzl=R \B) D =0, (66) 

where are the X-reflected versions of the VF-currents described above 

(i.e. are obtained from W s by repacing X — > —X in their expressions 

in terms of the fields (X,Y)). Since in general the ^-currents do not have 
symmetry with respect to these reflections, {W s } and {VFjP^}, as the sets 
of operators in ^"p, are different. These sets have some intersection though, 
as we will explain in the next subsection, and further exploit in Section f5. 21 

4.2 More on the W-algebra 

Here we want to make few statements on the structure of the VF-algebra of 
the hairpin model. 

The easiest way to generate the higher-spin VK-currents is to observe that 
the "screening operators" associated with the exponentials ()58(l commute 
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with the parafermionic currents 



9(z) = {V^+2 d z Y +iy/nd z X) qJ^ Yr{z) , (67) 
= ( Vn + 2 d z Y -i^ d z X)e~^^ YR(z) , 
in the same sense as the ^-currents do, i.e. 

j dw V(z) V±(w) = 0, jf du; 9*(z) V±(w) = 0. (68) 

In (|67[) the in the exponential stands for the holomorphic part of 

the local field Y(z,z) = Yr{z) + Yl(z), and therefore the fields (|6*7|) are 
not local - they extend the notion of the parafermions of |42| to non- 
integer k = —n — 2. Nonetheless, both ^ and ^* are local with respect 
to the exponentials (|58j) . hence the integration contour in (|68|) - a small 
contour around z - is indeed a closed one. It follows from (|68f) that all the 
fields generated by the operator product expansion of ^f(z)^>*(w) satisfy 
Eqs.(JS7J) 35 j. Thus we have 

(z)tf*(0) = z~—z { - — 1_ - '± (T(z) + r(0)) + (69) 

(W 3 (z) + W 3 (0))--^= {W i {z) + W i (0))+... }, 
\/n + 2 Ayn + 2 > 



where T and W3 are the same as in (j55j) and (|59|) . and the higher-order 
terms involve the higher-spin T^-currents. 

Explicit calculation shows that the current W4 appearing in (|69j) can be 
written as 



W A = iy 4 (sym) + d z 3 + aT 2 + b d 2 z T , (70) 



where 



w (sym) = (4n 2 + ^ + 4) (9 2 X)2 + ^2 + 7n + 2 ) (02y)2 + 

n (3n + 4) (d,X) 4 + (n + 2)(3n + 2) (d,Y) 4 + 

6n(n + 2) (9 z X) 2 (^y) 2 , (71) 

and O3 is some local field. The last two terms in (|70J) are irrelevant - they 
represent the ambiguity in adding composites of the lower VF-currents - and 
exact values of the constants a, b are not important here. Explicit form of 
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the field O3 is not important for the present discussion either. Important is 
the symmetry of W4 with respect to the X — ► —X reflection. 

Similar property can be observed for the next few VT-currents of even 
spins. Using a freedom in adding derivatives and composite fields built from 
the lower-spin VF-currents, the fields W 2 k can be brought to the form 

W 2k = W^ m) + d z 2k _ x , (72) 

where the expressions for m) are even with respect to the reflection 
X — > —X. We believe this statement is valid for all even-spin W-currents. 
Since the "right hairpin" ^-currents in (|66|) are X — > —X reflections 

of the W s , this property can be restated as follows: the even-spin VK-currents 
can be defined in such a way that 

^^i? m od Z (...), (73) 

where the notation wj c ^ = W s is used to stress that the currents W s de- 
scribed in this section generate the symmetries of the left hairpin. 



4.3 Short-distance expansion of the paperclip 

Qualitative discussion at the beginning of this section suggests that the 

hairpin overlap amplitude B(P,Q) controls the k — > asymptotic of the 

paperclip partition function Z{ P, Q \ k ). More precisely, B(P, Q) gives right 

asymptotic form of Z( P, Q \ k ) when P lays in the upper half-plane 9m P > 

0, and one has to take B(—P, Q) to describe the asymptotic of Z(P, Q | k) 

when P is in the lower half-plane. Since B(P, Q) vanishes in the limit R — > 

[ p 

if P is taken in the "wrong" half-plane (note the factor k ^ in ()65j0 . this 
in turn suggests that the overall R — * asymptotic of the partition function 
is correctly expressed by the sum 

Z(P,Q\k)\ k ^ ^ B(P, Q) + B(-P, Q) . (74) 

What can be said about corrections to this leading asymptotic? Let us 
assume that k is small but not asymptotically small; that means that the 
X-size of the paperclip in Fig.l is large but finite. Also assume that QmP 
is positive. Then the functional integral (|16|) is still dominated by the field 
configurations (X, Y) with X close to the left end of the paperclip, with 
some fluctuations towards its right end. One then expects to have two kinds 
of corrections. When the fluctuation is sufficiently small, the functional in- 
tegral "feels" only the small deviations of the shape of the paperclip @ from 
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the shape (|53|) of the left hairpin. This leads to the perturbative corrections 
to the above leading k — > asymptotic; simple dimensional analysis shows 
that these corrections behave as integer powers of k«, or equivalently, as 
integer powers of r 2 . This part of the structure is already explicit in the 
semiclassical expression (jSHJ- Also, there are large fluctuations, sensitive to 
the fact that the right end of the paperclip allows for the passage from one 
leg of the left hairpin to another. Important example of those are the instan- 
ton fluctuations discussed in Section l3~2l The instantons generate powers of 
r n or, equivalently, integer powers of k. 

This structure is neatly captured by the form 

Z(P,Q\k)=B(P,Q) F(P,Q\k)+B(-P,Q) F(-P,Q\k), (75) 

where the function F(P,Q \ k), apart from the overall factor k 2k (see Section 

2 

I2.2|) . is a double power series in integer powers of and n, 

oo 

F{P,Q\k) = k 2k fij(P, Q ) ^ " . (76) 

with _fo,o = 1- The exact splitting into two terms in (|75|) is not easy to justify 
on general grounds, but can be supported by the following arguments. First, 
recall that such splitting in the semiclassical expression (|36j) corresponds to 
isolating contributions from two saddle points, the "right" one near the left 
end of the paperclip, and the "wrong" one near its right end (We are still 
assuming that QmP > 0; otherwise the above qualifications reverse. They 
have nothing to do with the authors political standings.). Similar splitting 
occurs in the semiclassical one-instanton contribution (see Section [3.2(1 . for 
the same reason. More importantly, the full expression has to take care of 
the poles of the factor B(P, Q) - recall that Z(P,Q \ k) must be an entire 
function of P. The expression (|75j) is the simplest form fit for this job. The 
poles of B(P,Q) are located at the points in the lower half of the P-plane 
where i^/nP or i ^ take non-negative integer values. At this points the 
i— 

factor k in B(P, Q) "resonates" with certain terms of the expansion (|76|) 
in the second term in (|75j) . Therefore, the poles of the factor B(P,Q) in 
the first term can (and must) be cancelled by poles in appropriate terms of 
the expansion of the second term. For the "perturbative" poles at i y/nP = 
0, 1, 2 . . . this mechanism is evident in the semiclassical expression (|36|) . and 
it can be confirmed for the first of the "instanton" pole at i — ^ = 1 by the 
form of the one-instanton correction (|52|) . The form similar to (|75|) , together 
with this mechanism of the pole cancellation, was previously observed in the 
boundary sinh-Gordon model |43| . 
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5 Integrable Paperclip 



In this section we would like to argue that the paperclip model is integrable. 
We will put the arguments in turn below, but let us first formulate the prob- 
lem for our case, where the bulk theory is free CFT and all the interactions 
occur at the boundary. 

First, we need to have integrable bulk theory. This is not a problem 
since the bulk theory is free, as described by the action (^Q). The free theory 
certainly has an infinite set of commuting local integrals of motion. Actually, 
it has several such sets, as we are going to see below. The existence is clear, 
but to the best of our knowledge such sets have never been completely 
classified. The problem is relatively easy to address, and we believe it still 
can have surprises in store. 

For the sake of this discussion it is convenient to change to the cylindrical 
coordinates (v, v) = (r + ic, r — kr) related to the disk coordinates (z, z) 
through the logarithmic map 

z/R = e iv / R , z/R = e^ /R , (77) 

because relevant integrals of motion are going to be homogeneous in this 
frame. The disk \z\ < R becomes a semi-infinite cylinder r = r + 2irR, 
a > 0. The action still has the form (^Q), with the derivatives d z ,dz replaced 
by d v ,dy, and the integration limits changed accordingly. The fields X still 
obey the constraint (J2J) at the boundary, which is now placed at a = 0. The 
effect of the exponential insertion in © is accounted for by imposing an 
asymptotic condition on the behavior of the field X at the infinite end of 
the cylinder: X — > i ^ P as a — > +oo (iP real). As usual, the bulk equations 
of motion 9(9 X = suggest that <9X is a holomorphic field, and hence any 
composite field built from 9X and the higher derivatives is a holomorphic 
field as well 7 . We will generally denote such polynomial as P(v). Assuming 
that the coordinate a along the cylinder is interpreted as the Euclidean time, 
the integral l[P] = § ^ P(v) over closed contour around the cylinder gives 
rise to a conserved charge of the bulk theory. 

Now, let's assume that P s +\(v ) is an infinite sequence of the polynomials, 
such that all the associated charges commute, i.e. [I S ,I S /] = 0, where 
I s = I[P s _|_i] (we will come back to the question of what it takes to have such a 
sequence). Also, let P s+ \{v) be the corresponding sequence of "left-moving" 
currents, with 8 replacing d, and I s = I[P s+ i] = § |^ P s+ i(v) are the "left- 
moving" charges; the I s / also commute among themselves and with all I s . 

7 Here and below in this section d and 8 will stand for the derivatives over the cylindrical 
coordinates, d = d v and 8 = dv- 
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The theory with the boundary at a = is integrable if the boundary state 
| B } satisfies the equations Q for all members of the sequence. Roughly 
speaking, the meaning is that the boundary neither emits nor absorbs any 
amount of the combined charges I a — I s . The equations Q are the property 
of very special "integrable" boundary conditions, which have to be such that 
the differences P s +i(v) — P s +i(v), when specified to the boundary a = 0, 
reduce to total derivatives d T (. . .) • 

With this understanding, there are in principle two opposite routes of 
approach to the integrable boundary interactions. One is in direct analysis 
of a given boundary interaction, having in mind proving (or disproving) 
that the system of local integrals satisfying @ exists. Another is to start 
with the search for suitable system of commuting integrals of motion {I s }; 
once that is found, one can try to identify local boundary condition which 
is compatible with Eqs.@. In the next section we describe what is likely to 
be the commuting sequence {I s } associated with the paperclip model. 

5.1 Paperclip integrals 

The polynomials P(v) can be classified according to their spins, which coin- 
cide with their scale dimensions and equal to the total number of derivatives 
d involved. Below we assume that the subscript s + 1 indicates the spin of 
the field P s +i(v). We will look for a sequence of polynomials of even spins, 
{P2(v), Pa(v), . . .}, such that all the operators I s = § |^ P s+ \{v) commute 
among themselves. Having in mind the paperclip model, we assume that 
P2 = T, where T coincides with the energy-momentum tensor of the bulk 
CFT, i.e. 

h = -f^[(dX) 2 + (dY) 2 ], (78) 

so that the difference Ii — Ii coincides with the momentum in the T-direction. 
We also assume that all the higher P s respect obvious symmetries of the 
paperclip ((2j), i.e. they are symmetric with respect to both the reflections 
X — > —X and Y — > —Y. Then one can try generic polynomials with these 
symmetries for P^{v), Pq(v), . . . with the coefficients to be determined from 
the commutativity conditions [I s , I s /] = 0. The equations [Ii , I s ] = 
are satisfied identically for any polynomials P s +\(v) with constant (i.e. v- 
independent) coefficients. But the first nontrivial equation [13,15] = 
turns out to be surprisingly rigid. Besides having few "trivial" solutions 8 , it 

8 The "trivial" solutions are as follows: One is P 4 = T 2 , P 6 = T 3 - | (dT) 2 (all 
composites here are in terms of the Virasoro OPE); it corresponds to the "KdV series" 
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fixes uniquely a one-parameter set of solutions. Using suggestive notation n 
for the parameter, the corresponding integral I3 is written as 



6 J 2vr L 6(3n + 2) V ; (3n + 2)(3n + 4) V J y 1 

W^TT) {dY) + 6(3n + 2)(3n + 4) ^ + (79) 

±n 2 + 7n + 2 - ; 
6(3n + 2)(3n + 4) 1 ; 

(the expression for I5 is too long to be comfortably written down here and we 
present it in Appendix). The integral (|79|) is a special case of the first non- 
trivial integral of motion written down in |44j in relation with the sausage 
model [2^ in dual representation. Once I3 is fixed as in (|79|). the equations 
[I3 , I s ] = with s = 7, 9, 11 . . . can be studied spin by spin. We have 
done that for s = 7 and s = 9; in these cases the above commutativity 
condition fixes the operators I s uniquely up to the normalizations. We take 
this as a strong indication that an infinite sequence of commuting integrals 
{I s ; s = 1, 3, . . . , 2k — 1 . . .} exists, whose first representatives are the oper- 
ators (fT5)) , (|79*)) . Our statement here is that the expression (f79*|) is essentially 
unique, once one demands it to be a member of the commuting sequence 
{hk-l}- 

The commuting sequence {l2fc-i} whose first representatives are (|78|) and 
(|79(l is the unique candidate for the system of local integrals of the paperclip 
model; anticipating direct relationship, we will refer to it as the "paperclip 
sequence". Note that besides the obvious reflection symmetries, Eq . (|79jl 
(as well as the expressions for the higher integrals I5, I7, . . .) respects rather 
subtle formal symmetry of the paperclip equation (j2J), 

X^Y, n^-n-2. (80) 

Also, in the limit n — > 00 the integral (|79jl (and hence all the integrals l2k—i 
of this series) becomes invariant with respect to orthogonal rotations of the 
vector X, in accord with the "circular-brane" limit of the paperclip. 

in which all the currents P s are built from the c = 2 Virasoro generator T (see Ref. [13] 
for details). The other solutions correspond to the KdV series with c = 1 involving only 
one component of the field X, either X, or Y, or one of the combinations X ± Y. We call 
these solutions "trivial" because their existence was assured upfront. 
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5.2 Paperclip integrals from PU-currents 

Important thing to note about the local integrals I s is that the expressions 
for their densities P s +i do not involve any particular scale. On the other 
hand, as was discussed in Section 0] in the short-distance limit k — > the 
paperclip becomes a composition of the left and right hairpins. Then, if 
one assumes that the paperclip boundary state | B ) satisfies the equations 
(jU, the same must be true for the boundary states of both left and right 
hairpins, 

{Is -%) |£) c = 0, (81) 

and 

(I s -I s ) \B) D =0. (82) 

This is indeed the case, at least for the first few integrals: the paperclip 
integrals I s "reside" inside both the PU-algebras, W^ c ' and of the 

left and the right hairpins, and the equations l)81j) and l)82j) follow from the 
equations 1)60)) and (|66|) . The way it works 9 follows from the properties of 
the even-spin VF-currents mentioned in Section 14.21 Using the cylindrical 
coordinates (v,v), the PU-symmetry equations can be written as 10 

[wP(v)-WP(v)} v= ^ \B) C =0, (83) 
[WpHv)-Wp\v)} v=d \B) D =0. 

Let us take even s = 2k and assume that the currents Wik are chosen in 
such a way that Ea. H73|) holds. Then, integrating l)83|) around the cylinder 

one arrives at Eqs.(jHU), © with operators I 2fc -i = § % W^ m) (v). This 
argument implies that the densities P2k{v) of the paperclip integrals coincide 
with w!^ m \v) up to normalizations. For k = 2 this fact is evident when 
one compares ()71(l with (j79|) ; it can also be checked directly for the next few 
integrals. If one recalls the definition of the VF-algebra through Eq. (|57|) . it 
also suggests that the "paperclip" sequence of local integrals {Ibfc-i} can be 
defined as the system of local operators which commute with four "screening 

9 The arguments in this section follow closely the ideas implicit in |44|. 
10 The fact that 183H are equivalent to IfciUH and IbbH is obvious when one chooses 
W3, Wa, ... to be the conformal primaries. It follows that these two forms are equiv- 
alent with any choice of the ^-currents. 
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charges", i.e. 



J V 

/ 

J V 



dwW, 



dwW, 



r(sym) 
2k 



r(sym) 
2k 





(84) 



where V± = V± are the exponentials 1)58(1 and V± are their X — > — X 
reflections 11 . 

6 Infrared Paperclip 

The R ^ oo limit of the paperclip model is governed by an infrared fixed 
point of the boundary RG flow, whose nature is not exactly known upfront; 
The infrared domain R —* oo is beyond access of the weak-coupling methods, 
no matter how large n is. When the RG time t grows, the paperclip in Fig.l 
first shrinks in the X-direction, and then collapses to something of the size 
~ 1. At this point its curvature is no longer small, and the semiclassical 
intuition based on the geometry of (J2J) is not too reliable. Nonetheless, the 
simplest idea, which corresponds to a naive extrapolation of © , © , is that 
the paperclip ultimately shrinks to a point, i.e. the infrared fixed point is 
the Dirichlet boundary conditions for both X and Y: 



Assuming this simplest scenario, and combining it with the integrability of 
the paperclip model, one can make predictions about the large-K expansion 
of Z(P,Q\k). Indeed, under this assumption the boundary-state operator 
B(k), Eg. 1)11)1 . is expected to have the asymptotic k — > oo expansion (|14j) in 
terms of the paperclip integrals I s . Since the vacuum overlap amplitude @ 
coincides with the vacuum- vacuum matrix element of B(k), this implies the 

11 This in turn suggests a possibility of "dual" description of the paperclip model, in 
which the boundary constraint is replaced by the boundary potential 



with Y being the T-dual of Y . Simple-minded idea of the coefficients A £1£2 being just 
constants does not seem to work - the coefficients must involve some additional boundary 
degrees of freedom akin to the quantum group generators in |45|. We did not yet work 
out all details of this dual theory, but we hope to come back to it in future. 



X B = (0,0). 



(85) 
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oo expansion 



log Z( P, Q | « ) ~ log(^) _ C 2fc _! k 1 -™ I 2k -i{P, Q) (86) 

fc=i 

of the logarithm of the function (JJJ) in terms of the vacuum eigenvalues of 
the paperclip integrals, 

l2fc_i |P) = P^ 2fc J 2fe _i(P,Q) |P). (87) 

The term log (5^) = — log(\/2) in Ea. (j8T)|) is the boundary entropy of the 
Dirichlet boundary condition (|85|). Despite the fact that the coefficients 
C*2fc-i are not known a priori, Eq.(|86|) gives rather strong prediction be- 
cause all the dependence on (P, Q) in (|8l)|) comes through the eigenvalues 
l2fc_i(P, Q). These eigenvalues are polynomials in P 2 and Q 2 which can be 
computed explicitly when explicit form of the paperclip integrals ^k-i ls 
known. Thus, using (|78jl. (|7T?|> and (|127j) one finds 

P 2 O 2 1 

j r P o) = npA + n{n + 2) p2Q2 + (n + 2) g4 _ 

31 J 96(3n + 2) 16(3n + 2)(3n + 4) 96(3n + 4) 
n(2n + 3)P 2 (n + 2)(2n + l) Q 2 ^ 18n 2 + 36n + ll 



and 

/ 5 (P,Q) 



48(3ra + 2)(3n + 4) 48(3n + 2)(3n + 4) 1440(3ra + 2)(3n + 4) ' 
n 2 P 6 (n + 2) 2 Q 6 



320(5ra + 2)(5n + 4) 320(5n + 6)(5n + 8) 

3n 2 (n + 2) P 4 Q 2 3n(n + 2) 2 P 2 Q A 



+ TTT7Z l w . TT + 



64(5n + 2)(5n + 4)(5n + 6) 64(5n + 4)(5n + 6)(5n + 
n 2 (3n + 4) P 4 (n + 2) 2 (3n + 2) Q 4 



64(5ra + 2)(5n + 4)(5n + 6) 64(5n + 4)(5n + 6)(5ra + 

3n(n + 2)(5n 2 + lOn + 2) P 2 Q 2 
' 32(5n + 2)(5n + 4)(5n + 6)(5n + 8) 
250n 2 + 225n + 30) P 
(5n + 4)(5n + 6)(5n 4 
(n + 2)(76n 3 + 206n 2 + 137n + 28) Q 2 
320(5n + 2)(5n + 4)(5n + 6)(5n + 8) 
1420n 4 + 5680n 3 + 7385n 2 + 3410n + 564 
20160(5n + 2)(5ra + 4)(5n + 6)(5n + 8) 



n(76n 3 + 250n 2 + 225n + 30) P 2 
+ 320(5n + 2)(5n + 4)(5n + 6)(5n + 8) ^ ^ 

+ 
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7 Solvable Paperclip 



In this section we propose an exact expression for the boundary amplitude 
The expression is in terms of solutions of special ordinary second- 
order differential equation. Similar expressions are known in a number of 
integrable models of boundary interaction, beginning with the work of Dorey 
and Tateo ^S]- Our proposal extends this relation to the paperclip-brane 
model. In this case no proof is yet available, but we will show in this section 
that the proposed expression reproduces all the properties of the paperclip 
amplitude described above. 

7.1 Differential equation 

Consider the ordinary second order differential equation 
d 2 



P 2 T^-7-(^-^)7T^+- 2 (l+eT 



q,( x ) = o, (90) 



dx 2 y 1 + e x ^ 4J (1 + e*) 2 

where p, q are related to the components of P = (P, Q) in JBJ), 

p=l^P, q=\^/n + 2 Q, (91) 

and k is proportional to the radius of the disk, n = E t R. With this iden- 
tification in mind, below we always assume that k is real and positive. In 
the semiclassical case b> 1 the parameters p, q here are the same as p, q 
in and n here relates to the renormalized r = r(R~ 1 ) in (J2J) through 
Eq.(j2J) within the two- loop accuracy 12 . 

Let ^-(x) be the solution of ()9U|) which decays when x goes to — oo 
along the real axis, and ^+(x) be another solution of (|90)1 . the one which 
decays at large positive x. We fix normalizations of these two solutions as 
follows, 

-» e KX as x -» -co , (92) 

and 



exp| - (| + «) x- k j — ((l + z)t-l)| as x^+oo.(93) 



12 The higher-loop terms in © are scheme-dependent, but it is plausible that a scheme 
exists in which the curve (2) is perturbatively exact, provided J3J is modified as follows 

k = (n + 1) r n (1 — r 2 ) . 

With this modification, the symmetry transformation 1801 of the equation JSJ leaves k 
invariant . 
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Let 

W[V + , = 4- -T- ( 94 ) 

dx dx 

be the Wronskian of these two solutions. Then, our proposal for the function 
(|28|). normalized in accordance with (|30j) . is 

tT (2 K ) 2K e- 2K 



To make more clear the motivations behind this proposal let us discuss some 
properties of the solutions of the differential equation (|§U|) . 

7.2 Small k expansion 

The differential equation (|9U|) has the form of one-dimensional Schrodinger 
equation with the potential V{x) defined by the last three terms in the 
brackets in (|§U|) . In our analysis below, it will be often convenient to split 
the potential into two parts, 

V{x) =y_(x) + y+(x), (96) 

where 

y_(x) = -p 2 - (q 2 - I) - &X , (97) 

v i f 1 + e x ^ 4J (i + e x ) 2 v ; 

and 

V+(x) =k 2 (1 + eT . (98) 

When k goes to zero the potential in (|9~U)) develops a wide plateau at 

1 « x « I log (4,) , (99) 

where its value is close to —p 2 . In this domain each of the solutions ^+(x) 
and x I'_(x) is a combination of two plane waves, 

*±(x)=D±(p,q) e +[ P x + D±{-p,q) e" 1 ^ . (100) 

Obviously, the Wronskian (|94|) is written in terms of the coefficients in ()100|) 
as follows, 

W [¥+, *_] -> -2ip (Z>+(p, q) £>_ (-p, ? ) - £>+(-p, 9 ) D_ (p, <?)) . (101) 
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Then leads to the k 
(f7i|l. with 



limit of the partition function of the form 



B(P,Q) = -2ip 



D-(-p,q) 



(102) 



where P, Q related to p, q through Eas.(|9"T]), The amplitude Z)_ (p, q) is easy 
to determine since for x <C i log f-4-) the term V+(x), Ea, (|98|) . is negligible. 
In fact, to handle properly the asymptotic condition (|92|) . it is convenient to 
retain the constant part k 2 of V+(x). With V+(x) replaced by k 2 , Eq. (|90|) 
reduces to the Riemann differential equation, and one finds 



x 



) = e [px (l + e- x ) [p - K 



(103) 



2F1 



+ k + q - ip,^ + k - q - ip,l + 2k; 



l+c- 



) 



as the At — > approximation of exact ^_(x) in the domain x <C ^ log (-tt). 
Note that the way H103|) is written, it captures the correction ~ k, but 
mistreats the terms ~ k 2 and higher. We will come back to this expansion 
later. From (|TU3|) . 



D_(p, g) 



r(2ip) 



g + ip)r(i + g + ip) 



(104) 



On the other hand, when x S> 1 the VL (x) part of the potential becomes 
a constant, while V+(x) can be approximated by K 2 e nx . To understand 
the quality of this approximation one can make a change of the variable 
x = Xq + ~ y, where xq = ~ log (■£-); the equation (|9"U1) then takes the form 

lb lb \ Zi r\i / ^^^^^^ 

J 2 



dy 2 



nr 



+ e 2 ^ + <yy (y) 



* = 0, 



(105) 



where 



<5V(y) 



as k — > 
2k - 



0, and indeed 5V admits expansion in powers of k« , 
4p2 _ 4g ,2 + 1 



ne 22/ + 



_2» „/ L 
e n + (J(Kn 



Therefore for x>l 

V 7rn L 1— ,l ' 
Here ^(2) is the conventional Macdonald function, and 

£ = 7£ + ^(-5), 



(106) 



(107) 



(108) 
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with 7£ = 0.5772... being Euler's constant, and if){z) = ^logT(z). The 
overall normalization in Eq. (|107|) is chosen to ensure the asymptotic form 
Hence 

D+0 ,, ?)= yz e£ «(3^ r( _ 2i ,). (109) 

Thus, Ea. (|102|) reproduces the correct form H65j) of the hairpin boundary 
amplitude B(P,Q). 

When k is small but finite, corrections to (|100|) can be obtained using 
perturbation theory. In computing $ ~{x) the term V + (x) is treated as the 
perturbation, and corrections appear in the form of a power series in k. In 
fact, the leading correction ~ k is already contained in ()103|) . and higher 
powers of k are generated by further iterations. On the other hand, 
is computed from (|105|) with the term 5V(x) taken as the perturbation. 
It is evident from (|106|) that corrections to *$>+(x) obtained through this 

2 

perturbation theory have the form of a power series in k». As the result, 
the partition function computed according to (|9*5l) appears in the form ((75|) 
with F(P,Q\k) being a double power series (|76|) . The low-order terms can 
be easily obtained by direct perturbative calculations outlined above, for 
instance, 



f d i 

h,o = ^ log 



(110) 



and 



T{\+e-q-\ P )r{\ + e + q-i P ) 
ip) - HI + Q - IP) + 7s + H ~ f ) + 2 log(2) - 2 



hl _ _ gi±i^zi) (»)i- x {m) 

n 4p 2 - 4g 2 + 1 \ r(l - \ - 2i H) 



2(V + i) ^ r(I-2H£ 



7.3 Semiclassical domain k< 1 <Cn 

2_ 2_ 

When k is small but n is large, so that k« ~ 1, all powers of k« in the above 
expansion have to be collected. Let us assume that r defined through © is 
not too close to 1, so that n (1 — r 2 ) ~ n. This regime corresponds to the 
semiclassical domain of the paperclip model considered in Section 03 

First, let us consider the case when p and q in ()90j) are of the order of 1; 
it corresponds to the case of "light" insertion in (|16|) . with (P,Q) ~ A=- In 
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this regime the term V+(x) in the potential has the effect of a rigid wall at 
some point x r , i.e. to the left from this point, for x r — x 3> -, the potential 
V + (x) is negligible, but to the right from it grows very fast, so that for x > x r 
the solution Vl/+(x) is essentially zero. More precisely, when x is below but 
close to x r the solution *S?+(x) is approximated by a linear function, 

^ +{x) k, a r {x r - x) . (112) 

However, the position of the wall x r and the slope a r depend on n in a 
somewhat subtle way. Here it is useful to trade k for the "running coupling 
constant" r defined through the equation k = nr n (1 — r 2 ), which is identical 
to the perturbative RG flow equation of the paperclip, Eq.©. Then 

x r =log(!^)+0(i), (113) 

and 

W 252 ^ ( i+o( ^)- (u4) 

We will explain these relations shortly. Using them, the Wronskian (|94|) 
is easily determined. Taking, for instance, any point x close to the left of 
the wall where both (jl(J3|) and Q112|) apply, one finds it equal to a r times 
the expression (jl03|) (with k set to zero) evaluated at x = x r - this leads 
exactly to (|35j) . Moreover, the one-instanton contribution (|52[) is also repro- 
duced when one takes into account the first-order term of the k expansion 
of Eq.JMJ). 

The case of "heavy" insertion can be handled similarly. The potential 
V+(x) still can be treated as a wall at x r , in the sense that at x r — x ^> ~ 
its effect is negligible, but it dominates at x — x r 3> \. Close to the wall's 
position, i.e. at 

\x- x r \ < 1 , (115) 

the part V-(x) is well approximated by the constant 

V-(x) f» V-{x r ) = \ n(l-r 2 ) 2 P r 2 , (116) 

where P r is exactly the expression ()43(l . On the other hand, in the domain 
(|115|) the part V + {x) behaves as the exponential 

V+(x) wn 2 (l-r 2 ) 2 e 2y , y = \ n (1 - r 2 ) (x - x r ) . (117) 
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Therefore, in the domain (|115|) the solution ^ + is approximated by the 
Macdonald function, 

^ rf K iv {2&), (118) 

where f = ^jj= . The normalization factor in front of i<C is fixed by matching 
lfTT5|) to the asymptotic condition (H3J) • For v -> and -y > 1 (fTTH|) 
reduces to (j!12j) . The Wronskian (|94|) can be evaluated in the domain 
- <C x r — x <C 1, where both (|103j) and (|118j) are valid; the result is exactly 
Q36JI with the extra factor (|4"il) added. 



7.4 Large «; expansion 

The case of large k is expected to describe the infrared limit of the paperclip, 
see Section |SJ For the differential equation (J90|) it is the domain of validity of 
the WKB approximation. Applying the standard WKB iterational scheme 
|46| one finds for the Wronskian (|94|). 



log W = log(2/s) + / dx \k(V{x) - V (x)) + 

J-oo I ° K 



1 (V'(x)) 2 

+ ... k(119) 



P 3 (x) 



where V(x) = k 1 y'V^x). The subtraction term 

n 

V (x) = (l + e x )* (120) 

in the integrand derives from the asymptotic conditions (|93|) . With the po- 
tential V(x) as in (|9b|). this generates asymptotic expansion of the partition 
function ([25]). 

oo 

logZ(P,Q| K )~ -ilog(2)-^ K 1 - 2fc / 2fc _ 1 (P,g), (121) 

k=l 

where l2k-i(P, Q) are polynomials of the variables P 2 and Q 2 of the degree 
2k. Their highest-order terms follow from the first term in the integrand in 

CEH, 



* n fc -' (n + 2)' r((fc-l/2)n + Q k ^ , 

£jj «(Jfc-l)! r((jfe-l/2)n + ife + Z) 1 j 
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which is in agreement with the highest-order terms of the vacuum eigen- 
values ((SB)) , lf%9*)) . It is also straightforward to generate the full polynomials 
evaluating the integral (|119f) order by order in kT 2 . This calculation re- 
produces exactly the eigenvalues (|55|) and (JHHJ), as well as few higher-spin 
eigenvalues (we have verified exact agreement up to ig). Note that Ea. (|122|) 
fixes normalizations of the operators I^fc— lj with this normalization all the 
coefficients C2/C-1 in (C s in (EU) are equal to 1. 



7.5 Circular brane limit 

As was mentioned in Introduction, in the limit when n —* 00 with the 
parameter g~ l = n (1 — r 2 ) kept fixed (so that k in Q has finite limit, 

k — > g~ l e 2 f ), the paperclip curve (J2J becomes a circle (|15[). It is interesting 
to see how this limit is obtained in the differential equation ()9(J|) . 

When n goes to 00 and k remains fixed, the term V+(x) in the potential 
(|96|) becomes infinite at any finite x. Correct limit is obtained by making 
first the shift of x, 

x = y — log n , (123) 

so that V+(x) = k 2 (l + ~ e^)" - . Then, taking the limit n — > 00 brings the 
equation (|9T)|) to the form 



d 2 P 2 

3-0 - -r eV + ft2 exp(e y ) 
dy 2 4 



*(j/) = 0, (124) 



where P 2 = P 2 + Q 2 . Thus, our proposal applies directly to the circular 
brane model ()15|) . with W- / ['I' + , being the Wronskian of two solutions of 
the differential equation (|124|) which are fixed by the asymptotic conditions 

^^(y)^e Ky as y -» -00 , (125) 

and 

f e y /"5 es dz ] 

-»■ expj - — -Ky-nJ^ ~^ Z ~ 1 )) as +00 -(126) 

More details about the circular brane model, and more about the properties 
of the differential equation l|124|) . are presented in a separate paper |26j . 
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8 Appendix 

Here we present an explicit form for the density P§ of the local integral 
I 5 = <f%P 6 (v): 

Pe = ~5(5n + 2)(5n + 4) ^ 
(n + 2) 2 {dY)6 



5(5n + 6)(5n + , 
,2 



3n 2 (n + 2) 4 



(5ra + 2)(5n + 4)(5n + 6) 

3ra(n + 2) 2 
(5n + 4)(5n + 6)(5n + 8) 

3n(4n 2 + 7n + 2) 
(5n + 2)(5n + 4)(5n + 6) 
3(n + 2)(4n 2 + 9n + 4) /a v\2 f s^v\2 
(5n + 4)(5n + 6)(5n + 8) 

3n(n + 2) (6n 2 + 13n + 4) , 2 , . , 



(5n + 2)(5n + 4)(5n + 6)(5n + 8) 

,2 



(dV) 4 (dY) 2 

{dxf {dY) 4 
(dxf {d 2 x) 2 

(8Y) 2 {d 2 Yf (127) 

(d 2 x) 2 (dY) 2 



3n(n + 2)(6n 2 + lln + 2) 2 2 



(dx) 2 (a 2 v) 



2 



(5n + 2)(5n + 4)(5n + 6)(5n + 8) 

12n(n + 2)(7n 2 + 14n + 4) ^ 2 ^ y , 2y 



(5n + 2)(5n + 4)(5n + 6)(5n + 8) 

(3 3 X) 



131n 4 + 575n 3 + 860n 2 + 500n + 96 3 2 



20(5n + 2)(5n + 4)(5n + 6)(5n + 8) 
131n 4 + 473n 3 + 554n 2 + 232n + 32 3 2 
20(5n + 2)(5n + 4)(5n + 6)(5n + 8) ' ' 
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